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Abstract 

Background: Microcystis aeruginosa is a species of cyanobacteria commonly found in a number of countries and 
frequently related to animal poisoning episodes due to its capacity to produce the cyanotoxin known as 
microcystin. Despite vast literature on microcystin structures and their deleterious effects, little is known about its 
synthesis by cyanobacteria. Therefore, this study used proteomic tools to compare two M. aeruginosa strains, 
contrasting them for microcystin production. 

Results: 2-DE gels were performed and 30 differential protein spots were chosen. Among them, 1 1 protein spots 
were unique in the toxin producing strain and 8 in the non-toxin producing strain, and 14 protein spots were 
shown on both 2-DE gels but expressed differently in intensity. Around 57% of the tandem mass spectrometry 
identified proteins were related to energy metabolism, with these proteins being up-regulated in the toxin 
producing strain. 

Conclusions: These data suggest that the presence of higher quantities of metabolic enzymes could be related to 
microcystin metabolism in comparison to the non-toxin producing strain. Moreover, it was suggested that the 
production of microcystin could also be related to other proteins than those directly involved in its production, 
such as the enzymes involved in the Calvin cycle and glycolysis. 
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Background 

Microcystis aeruginosa,^ worldwide species of cyanobac- 
teria, is frequently related to toxic water blooms. These 
blooms have occurred in many countries, such as the 
Philippines [1], Australia [2], the United States of America 
[3], Canada [4] and Brazil [5], causing water contamin- 
ation, risk of toxicity for hemodialysis and fishing prohib- 
ition. Moreover, numerous cases of animal poisoning due 
to cyanotoxins have been documented worldwide 
(livestock stupor in Australia [6], fish deaths in the 
U.K. [7] and cattle deaths in Switzerland [8]). 
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Nevertheless, the most serious case of human poi- 
soning attributed to cyanotoxins in drinking water, 
called "Caruaru syndrome", occurred in Brazil, in 
1996, when 76 patients died after hemodialysis treat- 
ment [9,10]. The deaths were attributed to the cya- 
notoxins microcystin and cylindrospermopsin. 

M. aeruginosa is able to produce the hepatotoxin 
microcystins [11], which are members of a remarkable 
family of more than 90 cyclic peptides that inhibit he- 
patic serine/threonine protein phosphatases PP1 and 
PP2A of terrestrial mammals [12]. Furthermore, it is 
well-known that environmental factors - such as 
temperature, light intensity, nutrients, salinity, pH - and 
culture age can induce changes in toxicity and toxin 
concentration [13-15]. Despite the M. aeruginosa 
genome having already been completely sequenced, 
the cyanobacteria metabolic pathway of microcystin 
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production has not been completely elucidated [16]. Other- 
wise, the microcystin biosynthetic gene cluster has already 
been characterized in Microcystis [17], Planktothrix [18], 
and Anabaena [19]. In this field, a proteomic approach 
could help to shed some light on toxin production. Never- 
theless, only a few studies of cyanobacteria have been con- 
ducted using proteomic tools, being basically restricted to 
the Synechocystis species [20-24]. Furthermore, most of 
these studies focus on a specific compartment, such as the 
thylakoid membrane [23], outer membrane [21], plasma 
membrane [22], periplasm and cytoplasm [25,26]. Other 
cyanobacteria species, however, have also been proteomi- 
cally considered, such as Anabaena variabilis [27], 
Anabaena sp. PCC 7120 [28,29], Nostoc punctiforme 
[30,31], Nevertheless these studies have been restricted to 
structural proteomic descriptions. 

Otherwise, comparative proteomics have also been 
used to find differentially expressed proteins of cyano- 
bacteria in response to environmental stresses like sali- 
nity [32,33], light [34,35], temperature [36], and pH [37]. 
Some studies have highlighted cyanobacterial changes 
becoming adapted to the symbiosis state [38,39]. How- 
ever, the majority of studies have been done with Syne- 
chocystis sp. PCC 6803. Environmental conditions, like 
higher or lower salinity, acidity, variation of C0 2 and Fe 2+ 
concentration are directly related to cyanobacteria develop- 
ment in the natural environment. Therefore, it is important 
to perceive how cyanobacteria respond to these environ- 
mental changes to better understand the bloom formation 
processes. 

Considering the enormous worldwide importance of 
M. aeruginosa in relation to freshwater quality, especially 
for public supply in reservoirs, it is interesting to note 
the low-quantity of research in proteomics [40,41] with 
comparative proteomics related to microcystin biosyn- 
thesis. To this end, this work aims to compare one nat- 
urally toxin-producing strain of M. aeruginosa with 
another naturally non-toxin producing strain using 
proteomic tools in order to better understand toxin 
production. 

Results and discussion 

Growth curves and microcystin detection 

Aiming to compare two M. aeruginosa strains according 
their ability to produce cyanotoxins, both were growth in li- 
quid medium according description previously described in 
literature. Although the differences are slight, there is cor- 
respondence between curves shown regardless of the wave- 
length used (Additional file 1: Figure SI). Microcystins are 
continuously produced during the cyanobacterial life cycle, 
but large quantities are commonly found at the end of log 
phase. For this reason, protein collection occurs at this 
growth point (Additional file 1: Figure SI, arrows) in which 



intracellular microcystin concentration is typically higher as 
previously observed for M. aeruginosa PCC 7820 [42]. 

In order to confirm if the M. aeruginosa strain PCC 
7820 is really toxin producing and also if NIVA CYA 43 is 
really a non-toxin producing strain, the presence of micro- 
cystins in both cultures were analyzed in CI 8 column 
reversed-phase HPLC [Figure la]. The fraction indicated 
by the arrow was eluted only from the M. aeruginosa PCC 
7820 extract after 7 min 48 s, with 40.7% acetonitrile. This 
fraction was concentrated, resuspended in a saturated solu- 
tion of a-cyano-4-hydroxycinnamic acid matrix and ana- 
lyzed by MALDI-TOF/TOF. The MS spectra showed that 
the fraction contained at least three different variants of 
microcystins [Figure lb]. Microcystin-LR (m/z 995) was 
the variant with the most intensive ions. The other 
detected ions were identified as [Dha]microcystin-LR (m/z 
981), [D-Asp 3 , Dha 7 ]microcystin-RR, microcystin-LW and 
[D-Asp 3 , ADM Adda 5 ] microcystin-LR, according to MS 
data comparison in the literature [43-45]. 

The observation of more than one microcystin variant 
was not surprising, because several cyanobacteria may 
be able to produce more than one microcystin variant 
[14]. Besides, in 50 to 75% of cyanobacterial blooms, 
toxicity is associated with the simultaneous production 
of several cyanotoxins [46]. 

Proteome of M. aeruginosa 

SDS-PAGE (Additional file 2: Figure S2) was used to 
evaluate the protein extraction protocol quality, the re- 
producibility of the protein preparations, and to quantify 
the protein extracts for preparative 2-D gel analyses. For 
2-DE, each strain was cultured in triplicate. From each 
biological replication, proteins were extracted from the 
whole-cell lysate. The M. aeruginosa PCC 7820 sample 
showed well-resolved protein bands with molecular 
masses between 10 and 116 kDa. Similar data was 
obtained from the M. aeruginosa NIVA CYA 43 sample 
(Additional file 2: Figure S2). In order to better investi- 
gate differential protein abundance visualized by SDS- 
PAGE, 2-DE analyses were performed. Protein extraction 
followed the same protocol used for SDS-PAGE, with 
the addition of a cleaning-up step, as described in the 
methodology. About 220 protein spots became visible 
and were clearly resolved in M. aeruginosa PCC 
7820 2-DE [Figure 2a] and about 202 protein spots 
in M. aeruginosa NIVA CYA 43 2-DE [Figure 2b]. The 
majority of spots were observed between pH 4 and 7, as 
previously described for Anabaena variabilis [27], 
Synechococcus sp. PCC 7942 [47], and for Synechocystis 
sp. PCC 6803 [23]. This was the protein pattern 
expected for cyanobacteria samples, although most 
cyanobacterial proteomic studies have used strips with 
narrower pi ranges (4-7) and/or subproteomes from 
fractionation methods [20,22,23,35,37,39,40,48]. These 
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Figure 1 (a) M. aeruginosa sample analysis in HPLC - reverse phase CI 8 column - and detection in 238 nm, following a non-linear 
gradient from 35% to 95% of acetonitrile/TFA 0.1%. The fraction indicated by an arrow was evaluated by MALDI-TOF^OF. (b) Mass spectra 
obtained by MALDI-TOF^TOF from M. aeruginosa PCC 7820 extract eluted fraction in HPLC and respective microcystin variant identifications. 



studies have provided valuable information, but for Micro- 
cystis, we thought it more appropriate to expand the pi 
range and protein solubility (by adding Triton-X-100) to 
provide a more global profile, considering that this is the 
first report that compares one naturally toxin-producing 
strain with another naturally non-toxin-producing strain of 
M. aeruginosa. 



Identification and functional classification of the 
differentially expressed proteins 

BioNumerics, software version 4.6, was used to analyze 
gel images in order to detect and select spots of interest. 
The five gel images obtained for each strain (biological 
replicates) were compared. At the end of the analysis, 
the software had provided volume values for each spot 
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Figure 2 (A) Proteomic map of M. aeruginosa PCC 7820. Gel 12.5%, stained with Coomassie, and pi 3-11. (B) Proteomic map of M. 
aeruginosa NIVA CYA 43. Gel 12.5%, stained with Coomassie, and pi 3-1 1. Grey, white and black arrows indicate, respectively, up-regulated 
proteins; down-regulated proteins and proteins exclusive to that strain. Dotted arrow was the spot control. 
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detected. With this information, linear regression ana- 
lyses were done, and thus r 2 was calculated. The linear 
regression presented r 2 = 0.75 for M. aeruginosa PCC 
7820 quintuplicates and r 2 = 0.70 for M. aeruginosa 
NIVA CYA 43. Linear regression analyses of technical 
replicate gels commonly showed correlations higher than 
0.80 [49, 50]. It is important to highlight that in this study, 
we worked with biological replicates, and therefore the r 2 
value was slightly lower. Theoretical Mr and pi were 
obtained from MASCOT analysis and compared with the 
experimental values that were estimated directly from the 
gel images, using BioNumerics software. Moreover, toxin- 
producing cyanobacterial protein spots were considered 
differentially expressed when present/absent in all repli- 
cates: when at least two-fold changed (up or down) com- 
pared to control (non-toxin-producing) and also when the 
standard deviation bars did not intercross. There was a 
well-defined correlation between the theoretical and ex- 
perimental Mr values as shown in Table 1. 

After the BioNumerics software analyses, 30 proteins 
of interest were chosen for identification based on differ- 
ential expression. These proteins were identified by 
MASCOT and corresponded to 21 different proteins 
classes [Table 1]. Almost all identified proteins showed 
sequence correspondence to M. aeruginosa. The only 3 
proteins which did not correspond to M. aeruginosa 
were mlO, ml 8, and m28; all of which showed similar- 
ities to other cyanobacteria. One spot (ml7), identified 
as superoxide dismutase, was used as a control to in- 
crease MS data reliability. Protein identifications showed 
that 61% of differential proteins are related to energy 
metabolism, 18% corresponded to membrane proteins, 
6% are carboxysome and periplasmic proteins, 12% are 
cytoplasmic proteins and 3% are hypothetical proteins. 
To avoid double counting, rubisco was considered only 
in energetic metabolism proteins, although this enzyme 
is located inside carboxysome. Four of these 21 proteins 
were found in more than one spot on the gel. Among 
these are the enzymes Rubisco (spots mOl, m02, mil), 
FBA class II (spots m05, m21), PRK (spots m08, m09, 
m29), and porins (spots ml2, ml3, ml4, m23) [Table 1]. 

Energy metabolism proteins 

It is interesting to observe that most of the differentially 
expressed proteins identified are related to energetic 
metabolism, and these proteins were up-regulated in the 
toxin-producing strain. It is important to declare here 
that we did not estimate the amount of new protein 
synthesis, but the protein abundance. This increased 
protein amount can result from higher synthesis (up- 
regulation) and/or stabilization of the identified protein. 
Ribulose bisphosphate carboxylase (rubisco), glyceralde- 
hyde-3-phosphate dehydrogenase (GAPDH), triosepho- 
sphate isomerase, D-fructose 1,6-biphosphatase class II/ 



sedoheptulose 1,7-biphosphatase (FBPase class II/SBPase) 
- fructose- 1,6-biphosphate aldolase (FBA), phosphoribulo- 
kinase (PRK), ferredoxin-NDP + oxidoreductase (FNR) 
were all up-regulated in the toxin-producing strain 
in comparison with the non-toxin-producing strain 
[Figure 3a]. Rubisco, GAPDH, PRK and FNR were also 
observed in thylakoid membrane of Synechocystis sp. PCC 
6803 the first time that thylakoid membranes from this 
cyanobacteria were used in proteomic studies [23]. Ru- 
bisco protein appears as isoform in 3 different spots with 
pi variation in the toxin-producing strain. FBA class II 
protein also appears in 3 different spots, with pi variation 
in the toxin-producing strain. However, this protein 
appears in only two different spots in the non-toxin-pro- 
ducing strain. The same result was observed for the PRK 
protein, which appears in 3 different spots with pi vari- 
ation. Nevertheless, the toxin-producing strain presented 
the same protein at an extra spot, with a molecular mass 
below that of the other two. 

The proteins identified here are all involved in the Cal- 
vin cycle, glycolysis and respiration: processes respon- 
sible for providing energy to the cell. In the toxin- 
producing strain, rubisco (mOl, m02 e mil), FBPase 
class II/SBPase (m06, m07 e m22), FBA (m05 e m21) 
and PRK (m08, m09 e m29) presented relative intensities 
of expression more than twice those seen in the non- 
toxin-producing strain. GAPDH (m04) was ten times 
higher in the toxin-producing strain. FNR was identified 
only in the toxin-producing strain. Therefore, it is sug- 
gested that overexpression of proteins involved in energy 
production by toxin-producing PCC 7820 probably 
means a compulsory requirement to produce microcys- 
tin, and this production may be synonymous of higher 
energy consumption. This has already been noted in 
studies with Streptococcus suis [51], Candida glabrata 
[52] and Yersinia pestis [53], which showed that expres- 
sion modifications in metabolic proteins suggest higher 
virulence. The proteome showed an increase in expres- 
sion in metabolism-related proteins, like fructose 
bisphosphate aldolase, glycerol-3-phosphate dehydrogen- 
ase and triosephosphate isomerase. According to these 
authors, an increase in glycolytic metabolism suggests 
higher virulence. Besides, changes in central carbon 
metabolism have been associated with virulence attri- 
butes and the inactivation of virulence regulators in C. 
albicans [54]. Proteomic results with the bacteria Yersi- 
nia pestis virulence indicate that expression changes in 
metabolic proteins may have a more significant role in 
virulence than was previously appreciated [53]. In the 
present study, it also seems that proteins involved in the 
Calvin cycle, respiration, and glycolysis play important 
roles related to microcystin production, since the Calvin 
cycle is the pathway most activated in the toxin-produ- 
cing M. aeruginosa PCC 7820. 



Table 1 Differentially expressed proteins identified from cellular extracts of axenic cultures of M. aeruginosa after tryptic digestion and MALDI-TOF/TOF analysis 



Spot n. Sequence 



Accession no Species 



Protein ID 



Sequence Mowse Function 
coverage (%) score 



Theoretical Experimental 
pl/Mw (kDa) pl/Mw (kDa) 



m01 r^FQGPPHGIWERD 

GIIIHIHR 
m02 r^FQGPPHGIWERD 

A\^~MGFVDLMR 

KXFEGPPHGIQAERD 
m03 KIAAEFPDNFRL 

IAAEFPDNFR 

m04 KGILEYNDLPLVSSDYRG 

m05 RKPTGEILAISRI 

YAGENFLRI 

RKPTGEILAISR 

KYAGENFLRH 
m06 KNIADHVAVEAMRE 

SIEELWWMDRPR 

FFHGGAR 
m07 LQEMGITNPDR 

FVDWHLFDQPK 

m08 RKPDFSAYIDPQRQ 
KVIVIEGLHPLYDERV 
RKPDFSAYIDPQRQ 
KVIVIEGLHPLYDERV 

m09 RKPDFSAYIDPQRQ 
KVIVIEGLHPLYDERV 
RGH^DDVMAAINSRK 
RKPDFSAYIDPQRQ 
KVIVIEGLHPLYDERV 
VLSVLLGMTIGQIHR 

m10 TLPDLIHSAPR 



166367530 
166367530 

166363998 

166365244 
166365988 



159027237 



196257787 



M. aeruginosa 
M. aeruginosa 



Ribulose bisophosphate 
carboxylase, major subunit 

Ribulose bisophosphate 
carboxylase, major subunit 



M. aeruginosa Ferredoxin-NADP 
oxidoreductase 



M. aeruginosa Glyceraldehyde-3-phosphate 
dehydrogenase 

M. aeruginosa Fructose-1 ,6-bisphosphate 
aldolase class II (FBA) 



166365743 M. aeruginosa D-Frutose-1,6-bifosfatase 
classe II / sedoheptulose 
1,7-bisfosfatase (FBPase 
classe 2/SBPase) 

166365743 M. aeruginosa D-Frutose-1,6-bifosfatase 
classe II / sedoheptulose 
1,7-bisfosfatase (FBPase 
classe 2/SBPase) 

1 59027237 M aeruginosa Phosphoribulokinase (PRK) 



M. aeruginosa Phosphoribulokinase (PRK) 



Cyanothece sp. FAD-dependent pyridine 
nucleotide-disulphide 
oxidoreductase 



5.0 
9.0 

5.5 

5.3 

49 



10.2 



18.0 



27.3 



82 C0 2 fixation (Calvin cycle) 

and photorespiration 

92 C0 2 fixation (Calvin cycle) 

and photorespiration 



44 Oxidation and reduction 

(ficobilissome, thylakoid 
membrane) 

101 Glucose metabolic process 
(phosphorilation) 

74 Glycolysis 



78 Carbohydrate biosynthesis 

(Calvin cycle) 



59 Carbohydrate biosynthesis 

(Calvin cycle) 



6.8 / 53.20 6.4 / 53.86 

6.8 / 53.20 6.2 / 53.84 

6.1 / 44.81 6.9/35.61 

6.3 / 36.64 6.6 / 38.49 

6.0/39.15 5.6/43.01 



4.7/37.66 5.1 / 42.99 



4.37/37.12 4.95/43.12 



1 07 Carbohydrate metabolic process 5.45 / 38.04 5.00 / 40.8 



1 1 2 Carbohydrate metabolic process 5.45 / 38.03 4.86 / 40.77 



25 Oxidoreductase 



6.41 / 46.87 8.41 / 42.9 



m11 



r^FQGPPHGIWERD 



166367530 



M. aeruginosa 



10.6 



91 



6.76 / 53.20 6.57 / 54.67 
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(Continued) 



RGITMGFVDLMRE 

KNHGIHFRV 

TFQGPPHGIWER 
ml 2 RDNDVTIDDLFYRA 

KFADINFVYVRS 

RLQAGNFNNTFNQSGPTRT 

KFADINFVYVRS 

LQAGN FN NTFNQSGPTR 

QSQTPAHQSRANQR 
ml 3 RDNDVTIDDLFYRA 

LQAGN FN NTFNQSGPTR 

RDVSPTAWAYEALRS 

KFADINFVYVRS 
m 1 4 LQAGN FN NTFNQSGPTR 

KFADINFVYVRS 

RLQAGNFNNTFNQSGPTRT 
ml 5 RGWINDINB-QNTONYPIIADGDRK 

KVIALSVDSAESHRG 
ml 6 QYFGTOFU/NLR 

RQYFGETDETVNLRV 

ml 7 RPDYISDFLTK 
YLDYQNR 

ml 8 IVQASTGLEVLSDSILVQKLR 

FLIATPLLLALMR 
ml 9 RLITYGWAGDVTPI EEIGLVGVRE 

YLSPGELDR 

DMDYYLR 
m20 RVPAEIVFDQGLGDLFVCRV 

m21 KYAGENFLRH 
RKPTGEILAISRI 
RKPTGEILAISRI 



166365540 



166365540 



166365540 



ribulose bisophosphate 
carboxylase, major subunit 



M. aeruginosa Porin type major outer 
membrane protein 



166366392 M. aeruginosa Peroxiredoxin 



15.7 



M. aeruginosa Porin type major outer 10.4 
membrane protein 



M. aeruginosa Porin type major outer 
membrane protein 



19.4 



1 66363870 M. aeruginosa Triosephosphate isomerase 1 2.2 



166368140 M. aeruginosa Superoxide dismutase (SOD) 9.0 



220908731 Cyanothece ABC transport related 5.8 
sp. PCC 7425 

166363768 M. aeruginosa Allophycocyanin alpha subunit 25.5 



166364254 M. aeruginosa Carbonic anhydrase, beta-type 8.3 
(CA) 

166365988 M. aeruginosa Fructose-1 ,6-bisphosphate 10.0 
aldolase class II (FBA) 



106 



C0 2 fixation (Calvin cycle) 
and photorespiration 



Transport (integral 
membrane protein) 



4.14/60.14 4.78/50.13 



96 Transport (integral 

membrane protein) 



84 Transport (integral 

membrane protein) 



Cell redox homeostasis 



4.14/60.14 4.64/50.1 



4.14/60.14 4.48/50.43 



4.17/16.55 4.68/25.39 



62 Triosephosphate isomerase 4.14 / 25.21 4.85 / 27.01 

activity (glyconeogenesis 
and glycolysis) 

40 Catalyzes of superoxide to 5.83/21.72 4.6/21.50 

peroxide and molecular 
oxygen 



78 Stress response 



1 04 Photosynthesis 
(phycobilisome) 



130 Carbon utilization 
59 Glycolysis 



6.07/64.52 4.99/18.61 

4.21 / 17.19 3.79/18.09 

8.45 /26.10 7.50/21.93 

6.00/39.16 5.39/43.1 



Table 1 Differentially expressed proteins identified from cellular extracts of axenic cultures of M. aeruginosa after tryptic digestion and MALDI-TOF/TOF analysis 

(Continued) 



m22 


RSIEELWWMDRPRH 
FFHGGAR 
FVDWHLFDQPK 
RFVDWHLFDQPKY 


166365743 


M. aeruginosa 


D-fructose 1 ,6- bisphosphatase 
class 2/sedoheptulose 1,7- 
bisphosphatase (FBPase 
classe 2/SBPase) 


14.2 


67 


Carbohydrate biosynthesis 
(Calvin cycle) 


4.37 / 37.68 


4.95 / 45.86 


m23 


DVSPTAWAYEALR 
DVERALTAPHLTR 


166365540 


M. aeruginosa 


Porin type major outer 
membrane protein 


4.6 


45 


Transport (integral 
membrane protein) 


4.14/60.14 


3.94 / 57.46 


m24 


RQYFGTOOTNLRV 
RQYFG^DOTNLRV 


166363870 


M. aeruginosa 


Triosephosphate isomerase 


6.5 


62 


Triosephosphate isomerase 
activity (gluconeogenesis 
and glycolysis) 


4.14/25.21 


4.85 / 27.01 


m25 


RTFTEVAPQPAPEPSVSPIRG 
RGEPQLN PG DYVLG RV 


166367494 


M. aeruginosa 


Hypothetical protein 
MAE_47530 


30.8 


69 


Non function related 


5.32/13.01 


8.39/17.08 


m26 


RAPYDESEWYHLDLYENKG 

FQQCLVQTR 

ROTVMSQAEIHSR 


159029497 


M. aeruginosa 


Phycobilisome rod linker 
polypeptide 


2.1 


119 


Photosynthesis 
(phycobilisome) 


9.24/ 32.21 


5.51 / 29.88 


m27 


LVSLGLLK 


159028959 


M. aeruginosa 


Cobalt ABC transporter, inner 
membrane subunit CbiQ 


3.1 


72 


Cobalt trasnport and 
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Figure 3 Differentially expressed proteins obtained from 2D gels of M. aeruginosa PCC 7820 (white) and NIVA CYA 43 (grey) strains 
related to (a) energy metabolism; (b) membrane proteins (external and internal); (c) periplasmic proteins; (d) cytoplasmic proteins and 
(e) hypothetical proteins. Vertical bars correspond to standard deviation. The asterisks indicate protein expression that differed significantly 
(p<0.05) between the toxin producing (PCC 7820) and non-toxin producing (NIVA CYA 43) strains by the Student's t-test. 



Carboxysome and periplasmic proteins 

Carbonic anhydrase (CA) and a protein with high se- 
quence similarity to zinc ABC transporter of Thermo- 
toga maritima MSB8 were identified and both were up- 
regulated by the toxin-producing strain [Figure 3c and 
Table 1]. The up-regulation of CA may occurs at the 
same time as the expression of zinc ABC transporter at 
the toxin-producing strain, since zinc is a metal required 
for the active site of CA, where this metal has a central 
role in carbon fixation for photosynthesis [55]. 

These results are in accordance with observations about 
energy metabolism proteins. CA is an enzyme that 



participates in the C0 2 concentrate mechanism (CCM). In 
aquatic environments with pH values above 7, inorganic 
carbon is available preferably in the form of carbonic acid 
(HCO3) [56]. Therefore, aquatic photosynthetic microor- 
ganisms developed CCM to generate an intracellular pool 
of inorganic carbon. 

Data reported here shows that CA was expressed by 
the toxin-producing PCC 7820 [Figure 3c], 1.8 times 
stronger than the non-toxin-producing strain. Probably, 
the toxin-producing strain has an improvement of car- 
bon fixation. This improvement could be directly related 
to microcystin production, since the use of carbons and 
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energy was necessary in this process. Moreover, this re- 
sult could also be related to an increase of metabolic flux 
through the Calvin cycle, mentioned above. In both 
cases, this reflects a more efficient carbon fixation that 
leads to PHA accumulation. With a higher quantity of 
carbon, it may be possible to maintain an efficient 
microcystin production. 

Cytoplasmic proteins 

The identified protein in spot mlO presented a high 
similarity with an enzyme from the flavoprotein family: a 
FAD-dependent pyridine nucleotide-disulphide oxidore- 
ductase, present in Cyanothece sp. This enzyme acts in 
cellular redox homeostasis processes and oxidation and 
reduction, where its molecular function is to carry out 
the activities of oxidoreductase and the electron carrier. 
Co-factors iron (Fe) and sulfur (S) are necessary for its 
action [57]. Flavoproteins seem to play an important role 
in bloom situations since they are frequently found in 
cyanobacteria collected from dense blooms. In environ- 
ments with dense blooms, Fe and copper deficiency 
often occurs [57]. Flavoproteins have the capacity to 
substitute ferredoxins in some cyanobacteria when they 
lack Fe, because these two protein classes are similar in 
relation to charge set and redox properties [58]. This 
protein was expressed by both strains at similar statis- 
tical levels [Figure 3d] suggesting that both cyanobac- 
teria may naturally have this advantage. Therefore, their 
growth may be favored in these environments, but toxin 
production may not necessarily be favored. 

Peroxiredoxin (ml5), observed only in the M. aerugi- 
nosa PCC 7820 gel, presents physiological significance 
as anti-oxidative stress in cyanobacteria [59]. This en- 
zyme catalyzes the reduction of various hydroperoxides 
[60] and seems not to be related to microcystin produc- 
tion. Otherwise, the tryptic fragments of protein spot 
ml 7 were identified as superoxide dismutase, an enzyme 
that protects the cell against damage caused by super- 
oxide radicals [61]. 

The PHA-specific enzyme acetoacetyl-CoA reductase 
(m30) was exclusive to the toxin-producing strain. In 
Synechocystis sp. PCC 6803, this protein has the function 
of PHA (polyhydroxyalkanoate) accumulation in the cell 
and, together with ketothiolase-p PHA specific enzyme, 
forms the first complete PHA biosynthesis route known 
in cyanobacteria [62]. PHAs are biodegradable polye- 
sters, carbon and energy storage compounds, which are 
deposited as insoluble inclusions at cytoplasm [62]. The 
fact that this protein appeared only in the toxin-produ- 
cingstrain indicates that to maintain the state of micro- 
cystin producer, M. aeruginosa PCC 7820 must have 
access to energy reserves. It may be considered that, 
since the toxin-producing strain has higher expression 
levels of proteins related to energy metabolism in 



relation to the non-toxin-producing strain, there is a 
need to produce proteins required for intracellular accu- 
mulation of energy produced. The expression of this en- 
zyme only by the toxin-producing stain suggests that 
producing microcystin means that more energy is spent 
due to the carbon and energy accumulation need seen in 
toxin-producing M. aeruginosa PCC 7820. 

Membrane proteins 

In this work, 4 spots were identified as external membrane 
porins. The porin protein appears in 4 different spots on 
the non-toxin-producing strain gel, with 3 isoforms and 
another spot with pi and MW different from the other 
three. Surprisingly, spots ml2, ml3 and ml4 were about 
twice as highly expressed in the non-toxin-producing in 
relation to the toxin-producing strain and spot m23 was 
exclusive to the non-toxin-producing strain [Figure 3b]. 
This observation represents an inversion of what has been 
seen for the other protein classes, and no toxin relation has 
been drawn until now. However, higher expression of porin 
in the non-toxin-producing strain seems to be in accord- 
ance with higher expression of peroxiredoxin. In relation 
to phycobilisome proteins, we identified allophycocyanin 
alpha subunit (ml9) only in the toxin-producing strain, 
phycobilisome rod linker polypeptide (m26) and phyco- 
cyanin alpha phycocyanobilin lyase related protein (m32), 
of which the last three only appear in the non-toxin-produ- 
cing strain. 

Spot ml 8 appeared only in the toxin-producing 
strain [Figure 3b] and presented high sequence similar- 
ity with an ABC transporter of Cyanothece sp. PCC 
7425. ABC transporters are a family of membrane pro- 
teins that carry out the transport of molecules rather 
than ions, and are involved in several physiological 
processes [63]. 

Otherwise, the ABC transporter was considered as a 
stress-related protein [36, 64]. In Spirulina platensis, this 
protein was detected in all subcellular fractions up-regu- 
lated when this cyanobacterium was submitted to high- 
temperature [64] and to low-temperature stresses [36] 
and its transcriptional expression pattern was well corre- 
lated with the protein expression pattern. 

Then, Amnuaykanjanasin and Daub [65] recently 
observed that the fungal Cercospora nicotianae disrupted 
at the ABC transporter gene {atrl) shows a dramatic re- 
duction in the production of the toxin cercosporin. 
These data [65] indicate that qualitative and quantitative 
differential expression of proteins in the mycelium pro- 
vides evidence to partly explain the mechanism of 
pathogen virulence differentiation. Similarly, our results 
indicate that the toxin-producing M. aeruginosa PCC 
7820 presents higher synthesis of proteins involved 
in metabolism in relation to the non-toxin-producing 
M. aeruginosa NIVA CYA 43. 
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Conclusions 

It have been observed that microcystin production pro- 
vides a photosynthetic and ecological advantage - 
adjusting rubisco synthesis [66], participation in pro- 
cesses associated with light reactions [67] and predation 
inhibition [68,69]. Putative functions of microcystin 
should result in a competitive advantage for those strains 
that are able to produce it. Therefore, microcystin non- 
toxin-producing strains need to develop compensating 
mechanisms. According Meipner et al. [70], non-toxic 
strains have genes coding for proteins involved in toxin 
production (peptide synthetases). Nevertheless, those 
strains seem to produce fewer gene copies in relation to 
toxics trains, mostly differing by genes content. Data 
reported here show that toxic strains are capable of pro- 
ducing proteins involved in energy metabolism at a 
higher quantity than nontoxic strains. However, in com- 
petition experiments, it was shown that nontoxic strains 
are better competitors for light than toxic strains [71]. 
These authors also suggest that non-toxic strains may 
invest their resources in other cellular functions since 
they are unable to utilize energy to produce microcystin. 
Our results point in same direction. Nontoxic strains 
probably have compensating mechanisms, using energy 
not spent in microcystin production in multiple cellular 
functions, like higher porin production and also in the 
improvement of cell size. Data reported here suggest a 
dualistic response: either microcystin stimulates the pro- 
duction of metabolism-related proteins or else metabol- 
ism proteins are expressed in higher quantities that 
enable the cells microcystin synthesis. In summary, en- 
ergetic metabolism proteins (Calvin cycle, glycolysis and 
respiration) were observed, being expressed in larger 
amounts by the toxin-producing strain (M aeruginosa 
PCC 7820). Corroborating these data, the expressions of 
other proteins that depend on energy production, which 
includes proteins related to carbon and nitrogen accu- 
mulations forming a kind of strategic cellular reservoir 
for eventual nutrient starvation, were also obtained. In 
an overview, the energy improvement observed in toxin 
production may provide a continuity of microcystin pro- 
duction under different environmental conditions. Add- 
itionally, the same proteins involved in stress response 
also seem to be linked to microcystin production, which 
included the ABC transporters and as well as the perox- 
iredoxin forms. This information leads to suggest the 
idea that an improvement of carbon and nitrogen fix- 
ation is necessary in the first moment. This modification 
may lead, in the next moment, to a higher ATP synthesis 
[Figure 4]. In parallel, the necessity of a stress response 
signal was also observed, correlated with the enhanced 
expression of cytoplasmic and membrane proteins. 
These processes together seem to orchestrate part of the 
microcystin synthesis process. Moreover, it is very 



interesting to notice that the Microcystis -cell machinery 
may work in a complex and harmonious manner, having 
the modulation of several proteins involved in microcys- 
tin production. Despite findings in energy metabolism 
and stress response provided here, it is also important to 
be clear that other processes could also be involved in 
cyanotoxin production that were not detected here by 
proteomic techniques. Another important issue is that 
we do not know at the moment if the non-toxin produ- 
cer has the microcystin synthetase (mcy) genes that en- 
code specific peptide synthetases, which have been 
shown to be involved in the construction of the toxic 
peptides. Nevertheless, studies performed by Cristiansen 
et al. [72] showed that gene cluster deletion is a rare 
event in nontoxic cyanobacteria. These modifications 
associated with the evolutionary diversification of the 
strain lineages lead to a combination of cyanobacteria 
that either lost or still contain the mcy-gene cluster in 
the same environment. Overall, proteomics technology 
may be considered a powerful tool used for uncovering 
the nature of microcystin production and its regulation 
in function of environmental changes. This understand- 
ing may significantly contribute to the development of 
compounds which may be employed directly in the en- 
vironment, aiming to block microcystin production by a 
toxin-pathway-specific means, thus preventing its release 
into water. 

Methods 

Organisms and culture conditions 

M. aeruginosa PCC 7820 and M. aeruginosa NIVA CYA 
43 were acquired from CSIRO (Commonwealth Scien- 
tific and Industrial Research Organization), Australia, 
and grown in a modified MLA medium (MgS0 4 .7H 2 0 
200 mM, NaN0 3 1 M, K 2 HP0 4 40 mM, H 3 B0 3 40 mM, 
NaHC0 3 200 mM, CaCl 2 .2H 2 0 200 mM, vitamins and 
micronutrients, pH = 7.5) in 1000 mL Erlenmeyer flasks 
(800 mL liquid). The cultures were maintained at 22°C 
under a 12:12 h light/dark cycle and illuminated with 
6 (imol of photons.m^.s" 1 of fluorescent white light. 
Cells were harvested from the growth medium at late 
exponential phase (OD 680nm = 0.7 - for M. aeruginosa 
PCC 7820, and OD 680nm = 0.8 - for M. aeruginosa NIVA 
CYA 43). The growth profile of M. aeruginosa was con- 
ducted in order to determine the end of the exponential 
phase. The growth curve was determined in 25 days, 
and the OD evaluation was conducted simultaneously 
every three days during this period. All measurements 
were done in triplicate at three different ODs: 680 nm, 
730 nm and 750 nm for each sample, considering wave 
length absorbed by cyanobacterial pigments, such as 
chlorophyll and bilins. These three different ODs were 
used to verify the correspondence between growth 
curves obtained by these different wave lengths and if 
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M. aeruginosa toxin producer strain PCC 7820 




M. aeruginosa toxin non-producer strain NIVA CYA 43 




Figure 4 Schematic representation of cellular differences observed in this report at proteomical level, between M. aeruginosa PCC 
7820 (A) and NIVA CYA 43 strains (C). Light microscopy images of M. aeruginosa PCC 7820 (B) and NIVA CYA 43 strains (D) were also provided. 
Blue arrows represent up- or down regulation processes observed. Yellow forms indicated the membrane proteins identified in this report. 



both strains grew equally well For cell extraction, 
OD 680nm was used to track cell growth. Although there 
were slight differences, correspondence between curves 
was shown regardless of the wavelength used (Additional 
file 1: Figure SI). For protein extraction, cultures of M. 
aeruginosa PCC 7820 and NIVA CYA 43 were centri- 
fuged at OD 680 nm = 0.7 and OD 680nm = 0.8, respectively, 
corresponding to the end of log phase. 

Microcystin detection 

During log phase, 1 mL aliquots of each strain were 
collected for evaluation. These aliquots were ground under 
liquid nitrogen (five times) followed by pelleting of insol- 
uble debris by centrifugation at 12,100 g for 5 min. Protein 
concentration was determined by Qubit™ Fluorometer 
(Invitrogen, Oregon, USA). Each aliquot was separately 
analyzed on a reversed-phase C-18 HPLC (GE Healthcare) 



column (Vydac 218 PT 5 \xm pore), at flow of 1 mL.min" 1 . 
This was run for 35 min and followed by a non-linear 
acetonitrile gradient (35% to 95%), using TFA 0.1% as ion 
pairing, according to Harada et al. [73] with modifications. 
The absorbance was monitored at 238 nm. The collected 
fractions were dried in speed vacuum (Labconco). Fraction 
analyses were conducted in a MALDI-TOF/TOF Ultraflex 
III (Bruker Daltonics). The chromatography fractions were 
resuspended in 20 \iL TFA 0.1% and mixed in a 1:3 v/v pro- 
portion with a saturated solution of a-cyano-4-hydroxycin- 
namic acid matrix. From this mixture, 1 \iL was manually 
spotted onto a MALDI target plate. The spectra were 
acquired in positive reflector mode. Ions with enough in- 
tensity were selected for cleavage by MS/MS. Microcystin 
fragmentation standard was generated by LIFT mode with 
external calibration using Peptide Calibration Standard II 
(Bruker). 
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Protein extraction and two-dimensional gel electrophoresis 

Cells were harvested by centrifugation (3,000 g for 30 min 
at 4°C) and mixed in buffer extraction (5 mM Tris-HCl 
containing 5 mM CaCl 2 .2H 2 0, 0.01% triton-X-100 and 
1 mM proteinase inhibitor cocktail (GE Healthcare)), fol- 
lowed by grinding under liquid nitrogen (five times) and 
pelleting of the insoluble debris by centrifugation at 
3,000 g for 30 min at 4°C. The protein concentration was 
determined using Qubit™ Fluorometer (Invitrogen, Ore- 
gon, USA) according to the manufacturers instructions. 
2-DE was performed with 13 cm gels with a 3-11 pH 
range loaded with 600 ug of protein. Prior to electrophor- 
esis, protein samples (biological triplicate for each strain) 
were TCA precipitated Pellets were washed three times 
with 100% cold acetone and cleaned with 2-DE cleaning- 
up kit (GE Healthcare), as indicated by the manufacturer. 
Precipitate was solubilized in 250 uL of an electrofocusing 
solution containing 8 M urea, 2 M thiourea, 2% (v/v) 
CHAPS, 1% (v/v) IPG buffer 3-10, 65 mM DTTand traces 
of bromophenol blue. Strips were rehydrated in this solu- 
tion overnight at 20°C. The isoelectric focusing of the pro- 
teins was performed on an IPGphor (GE Healthcare), and 
the running conditions were (i) 500 V for 1 h; (ii) 1000 V 
for 2 h; (iii) 3500 V for 5h30m and (iv) 3500 V for 10 h 
until a total of 40,000 Vh. The gel strips were then incu- 
bated in equilibrated buffer (Tris-HCl 1.5 M, 6 M urea, 
30% (v/v) glycerol, 2% SDS, 1% DTT, 1% iodoacetamide 
and traces of bromophenol blue) for 30 min before being 
pointed at the top of a 12.5% polyacrylamide gel and sealed 
with 0.5% agarose. Electrophoresis was carried out at 15°C 
and 67 mA/gel for 8 h. Proteins were detected by Coomas- 
sie Blue G-250, scanned using an image scanner (HP scan- 
ner Scanjet 8290) and evaluated with the BioNumerics 
software 4.6 version (Applied Maths NV, Belgium) as 
described below. 

Image analysis 

A total of five gel replicates for each strain (M. aerugi- 
nosa PCC 7820 and M. aeruginosa NIVA CYA 43) were 
selected and submitted to in silico analysis by BioNu- 
merics 5.1 software (Applied Maths). 16-bit TIFF gel 
image replicates (gray scale and 600 dpi) were aligned 
and screened by the software in order to detect the 
protein spots and determinate their concentration, mo- 
lecular weight (Mw) and isoelectric point (pi). After- 
ward, gel images of all strains were intercrossed and 
screened for spot differences and similarities. Statistic 
correlation analysis (R 2 ) was applied to evaluate differen- 
tial gel data. Correlations cut offs were applied and gels 
with R 2 lower than 0.7 were discarded. Mean value of 
spot volume (dpi) from spot densities was used for fur- 
ther comparison of protein expression, taking into ac- 
count spots with relative volume equal to or bigger than 
0.1 dpi. Moreover, the toxin-producing cyanobacterial 



protein spots were considered differentially expressed 
when present/absent in all replicates; when at least two- 
fold changed (up or down) compared to control (non- 
toxin-producing strain); when the standard deviation 
bars did not intercross. Finally, all spots that achieved 
the minimum requires were submitted to Students £-test 
and differences with p < 0.05 were considered significant. 

In-gel digestion and MALDI-TOF analysis 

In-gel digestion, including reduction with DTT and al- 
kylation with iodoacetamide, was performed according 
to Schevchenko et al. [74], with minor modifications, 
using 600 ng of trypsin gold-mass spectrometry grade 
(Promega) per spot and digested for 24 h at 37°C. The 
samples were air dried, resuspended in 10 uL of milli-Q 
water and analyzed in a MALDI-TOF/TOF Ultraflex III 
(Bruker Daltonics). One uL of each digest was mixed 
with 3 uL of a saturated solution of a-cyano-4- 
hydroxycinnamic acid matrix and was manually spotted 
onto a MALDI target plate. The spectra were acquired in 
positive reflector mode with 200 random shoots. Ions with 
enough intensity were selected for cleavage by MS/MS. 
Fragmentation standard was generated by LIFT mode 
with external calibration, performed using Peptide Cali- 
bration Standard II (Bruker) with accuracy of 0.3 Da. 

Database search 

Proteins were identified by using the MASCOT search 
engine (Matrix Science, London, UK) using the NCBInr 
database. The search parameters included carbamido- 
methylation of cysteines, oxidation of methionines, one 
miscleavage by trypsin and 80 ppm mass accuracy. The 
identification of the proteins was repeated at least once 
using spots from different gels. Identifications with prob- 
ability based on Mowse Scores > 50 were considered sig- 
nificant (P<0.05), and identifications with Mowse 
Score < 50 were then disregarded. 

Additional files 



Abbreviations 

ABC: ATP-binding cassette; Adda: (2S,3S,8S,9S)-3-amino-9-methoxy-2,6,8- 
trimethyl-10-phenyldeca-4,6-dienoic acid; CA: Carbonic anhydrase; 
CHAPS: [(3-Cholamidopropyl)dimethylammonio]-1 -Propanesulfonic Acid; 
OD: Optic density; D^: Dithiotreitol; FBA: Fructose-1,6-biphosphate aldolase; 
FBPase class ll/SBPase: D-fructose-1,6-biphosphate classll/sedoheptulose 1,7- 
biphosphate; FNR: Feredoxin-NADP + oxidoreductase; 
GAPDH: Glyceraldehyde-3-phosphate dehydrogenase; HC0 3 : Carbonic acid; 



Additional file 1: Figure SI. Growth curves of M. aeruginosa PCC 7820 
(a) and M. aeruginosa NIVA CYA 43 (b). Cyanobacteria development was 
evaluated at wavelength of 680 nm, 730 nm and 750 nm. Black arrows 
correspond to the end of log phase in which protein extraction occurs. 

Additional file 2: Figure S2. SDS-PAGE 15% by using crude protein 
cytosolic extract (50 ug.mF 1 ) from M. aeruginosa samples, (a) corresponds 
to M aeruginosa PCC 7820 and (b) corresponds to M. aeruginosa NIVA 
CYA 43. Gel was Coomassie stained. M corresponds to molecular mass. 
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HPLC: High performance liquid chromatography; MALDI-TOF: Matrix assisted 
laser desorption ionization time of flight; NRPS: Non ribosomal peptide 
synthesase; PHA: Polyhydroxyalcanoate; pi: Isoelectric point; PKS: Polyketide 
synthase; PRK: Phosphoribulokinase; 2-DE: Bidimensional electrophoresis; 
DSD-PAGE: Sodium dodecil sulphate-poliacrylamide gel electrophoresis; 
SOD: Superoxide dismutase; Rubisco: Ribulose bisophosphate carboxylase; 
TCA: Trichloroacetic acid; Vh: Volt hours. 
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